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ABSTRACT 

A numerical procedure is presented to calculate transmittance degradation caused by contaminant films on spacecraft 
surfaces produced through the interaction of orbital atomic oxygen (AO) with volatile silicones and hydrocarbons from 
spacecraft components. In the model, contaminant accretion is dependent on the adsorption of species, depletion reactions 
due to gas-surface collisions, desorption, and surface reactions between AO and silicone producing SiO x (where x is near 2). 
A detailed description of the procedure used to calculate the constituents of the contaminant layer is presented, including the 
equations that govern the evolution of fractional coverage by specie type. As an illustrative example of film growth, 
calculation results using a prototype code that calculates the evolution of surface coverage by specie type is presented and 
discussed. An example of the transmittance degradation caused by surface interaction of AO with deposited contaminant is 
presented for the case of exponentially decaying contaminant flux. These examples are performed using hypothetical values 
for the process parameters. 


1. INTRODUCTION 

Surface contamination of components continues to be a primary concern with regard to maintaining performance levels of 
spacecraft in low Earth Orbit (LEO). 1 * 4 As recent as the 8 lh International Symposium On “Materials in a Space Environment” 
held in Arcachon, France during June 2000, it was determined to be the problem of most widespread interest according to a 
survey response. Particular contaminants of interest include silicones and hydrocarbons released from spacecraft 
components. 5 Silicones are widely used as potting compounds, adhesives, seals, gaskets, hydrophobic surfaces, and atomic 
oxygen (AO) protective coatings. Some of these silicones are volatile. Volatile silicone molecules can arrive at and deposit 
on exposed spacecraft surfaces along with volatile hydrocarbons. Some of these surfaces are also exposed to AO that 
interacts with and oxidizes deposited silicones and hydrocarbons. It is believed that some of the oxidation products contribute 
to the formation of contaminant films, while others are volatile products. Contaminant films pose a risk to spacecraft 
performance because they can act as an optically absorbing film on the surfaces of sun sensors, star trackers, and optical 
components or increase the solar absorptance of thermal control surfaces. This paper addresses the attenuation of optical 
sensor transmittance caused by an accumulation of surface contaminants formed by the reaction of in-space AO with volatile 
silicones and hydrocarbons discharged from spacecraft materials. 

AO is abundant in the LEO environment and readily reacts chemically with many exposed organic materials during a 
mission. A chief product of such a reaction with silicone-containing materials is SiO x . To help predict the effects of the space 
environment interaction with organic compounds to form contaminant film on optical surfaces, a numerical model was 
constructed. In the model, specified silicone and hydrocarbon molecular fluxes strike a fused silica substrate and are partially 
adsorbed. Some of the adsorbed species chemically react with a concurrent specified AO flux. The model is based on a 
double-layer (surface layer and first sub-layer) bookkeeping procedure that relates the respective fractional surface areas of 
the participating species to parameters including volatile species flux, surface adsorption and desorption coefficients, and AO 
chemical reaction probability. Because of sublation through desorption and other depletion mechanisms that partially expose 
the sub-layer, the fractional species coverage is also book-kept for the first sub-layer of adsorbates. Updating the relative 
surface coverage is dependent on the fractional surface and sub-layer compositions, on current species fluxes, and on the 
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probabilities of adhesion, desorption, depletion and chemical reaction. Because surface adsorption even with only one species 
type is a complex process, and because the detailed interaction mechanisms are not clearly understood or readily calculated, 
many processes are treated in a simplified manner and include the use of hypothetical values for the process parameters. 

Surface reactions are attributed to AO and silicone molecules present on the surface. Reactions can also occur due to the AO 
flux and silicone flux impinging on the surface sites occupied by these species types. It is assumed that the molecular SiO x 
product forms a surface coating which, depending on species coverage, traps other molecules beneath it. The trapped 
molecules that may influence optical absorptance are volatile silicone or hydrocarbon molecules. The contaminant film 
growth depends on the surface reaction of adsorbed AO and silicone molecules, on the reaction of AO flux with adsorbed 
silicone molecules, and on the reaction of silicone flux with adsorbed AO atoms. The accumulated contaminant layer is 
computed by modeling the accretion. The variation in optical transmittance (neglecting UV darkening effects) can be 
calculated at successive stages by comparing the surface portions of contaminant film due to accrued SiO x versus the surface 
portion due to trapped hydrocarbons and silicones. 

2. ADSORPTION AND DESORPTION CONCEPTS 

Depending on the energy of adsorption, surface adsorbate bonding is classified under two broad types: physisorbed or 
chemisorbed. 6 ' Physisorption is associated with weaker intermolecular forces such as van der Waals interactions, whereas 
chemisorption is due to stronger covalent, ionic and metallic bonding caused by electron exchange between adsorbate and 
adsorbent. The weak physisorption bonds allow the adsorption process to be easily reversed, and hence average residence 
times for physisorbed molecules is short. Because adsorption is an exothermic process, the residence time, T r , of a 
physisorbed molecule depends on the amount of excess energy it retains after being adsorbed. The residence time for a 
physisorbed or chemisorbed molecule can be related to the typical vibration period x 0 (-10' 13 s) and to adsorption enthalpy, 
A h a , by the following Arrhenius type expression 


T,. = Te 


-AhJRT 


[ 1 ] 


where R is the gas constant and T is the substrate temperature. Upon examining this expression, it is clear that the lifetime of 
an adsorbed molecule depends strongly on the potential well into which it is adsorbed and on the temperature. For example, 
at a temperature of 273 K and using energies of adsorption equal to -30 kJ/mole and -lOOkJ/mole for physisorption and 
chemisorption, respectively, the corresponding residence times of 5.5xlO" 8 s and 1.4xl0 6 s are obtained. These times are 
many orders of magnitude apart. Furthermore, increasing the temperature to 373 K correspondingly results in decreasing the 
residence times to 1.6xl0~ 9 s and 10 s, respectively. Often for a process, it can be assumed that equilibrium is achieved 
between adsorption and desorption over the time period of the chemisorption process. For our purposes, the chemisorption of 
silicones can be considered as being carried out in successive steps, as physisorption followed by activation to a bound state 
through its interaction with either adsorbed AO or AO flux. Desorption is also an activated process; here the flux of any 
component may activate desorption. In this paper the microscopic details of the mechanisms are not evaluated but 
macroscopic procedures are used based on probabilistic arguments to represent the gas-surface processes. 

Because the typical time scale of physisorption is much smaller than the time of contamination (months or years in LEO) the 
details of individual atoms or molecules striking a surface are replaced by a net flux of molecules that, based on 
accompanying sticking probabilities, produces a coverage representative of a detailed balance approach. Consequently, the 
small time increments associated with microscopic events are linked together into much larger steps, which in turn may be 
small compared to the total experimental period. The sticking probability, P 5 , for an adsorbate can be defined descriptively as 

p _ rate of adsorption of molecules at the surface ^ 

rate of collisions of molecules with the surface 


It is understood that P s itself may depend on coverage since it is dependent on A h a , which depends on surface coverage. 
Thus, the rates of adsorption and desorption are highly influenced by the activation energy barrier, which depends on surface 
coverage. At a particular instant, the existing surface will be covered with various species that may be randomly distributed 
or perhaps distributed as clusters of various single specie islands. Although the particular arrangement of adsorbed molecules 
affects the thermodynamics of the process, effects of site arrangement are not addressed in this paper. 
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3. FRACTIONAL SURFACE COVERAGE AND SUB-LA YER COVERAGE 


For conceptual purposes it is convenient to group the various surface species into contiguous regions. This is also useful in 
bookkeeping the different processes. Accordingly, a schematic drawing of the surface coverage two layers deep is given in 
figure 1. This schematic indicates that the upper surface is composed of the four species, with the fractional coverage given 
by F g , F s , F hc and F ao for the silica-glass, silicone, hydrocarbon and AO, respectively. These are the four larger rectangles that 
constitute the unit area. The coverage components of the first sub-layer are also indicated and denoted by F y, where both 
indices run over the species range. Thus, for example, the second small rectangle from the left along the top row is denoted 
F s s and indicates the fraction of glass in the top layer that is covering silicone. From this definition, it follows that the sum of 
F 8 g , F 8 5 , F 8 hc and F s ao equals one because they constitute the four sub-layer fractions covered by glass. These fractional 
components are multiplied by F g to obtain equivalent sub-layer fractions of the total unit area (not just the glass area). This 
holds for the fractional coverage below the other species. The connection between the sub-layer fraction in terms of coverage 
specie and the sub-layer fraction in terms of total area is given by the following equation 

f;=F,F‘ with X F 1 =1 and £ f)=F i> [3] 

j j 

where /y has been introduced to denote the fraction of the total unit area that sub-layer type i covered by type j represents. In 
what follows, the F, and f j are referred to as absolute fractions because they are fractions of the total unit area. Similarly, the 
F l j are called relative fractions because they are fractions of the non-unit F,. 

4. RELATIVE CHANGE OF COVERAGE BASED ON THE PROCESS 


4.1 Adsorption 

The fractional area of coverage can be computationally modeled by approaching it from the perspective of exposing a surface 
to a fluence burst of m molecules per unit area. The following processes are based on a unit time period. In a subsequent 
section, the fractional coverage will be expressed in iterated form. Based on a sticking probability, P s , and n available 
discretely defined sites per unit area, the fractional coverage, 6 , is derived from simple statistical considerations, when 
allowing for multiple hits on a site, as 


e=\- e ~ {m ^ /n) . [4] 

This provides a convenient form from which to obtain the change in species coverage given m , n, and P s . For given surface 
conditions (fixed n and P s ), it can be seen from Equation [4] that 6 approaches unity for large atomic or molecular fluences. 
Figure 2 shows curves of fractional coverage versus fluence-to-site-density ratio for various P s . At higher values of P v , the 
curves become more nonlinear, indicative of a larger percentage of molecules hitting occupied sites. Equation [4] will be 
used to determine fractional coverage based on adsorption. It is assumed that P s is known or prescribed for every possible 
species-flux-on-surface combination. 

Based on the preceding discussion, the relative fractional coverage due to the adsorption process is defined for the various 
flux-surface combinations. The relative fractional coverage due to each of the other processes is discussed subsequently. In 
this model, the three species of fluxes are AO , silicone, and hydrocarbon, which can land on sites occupied by the other 
species as well as on the silica glass substrate. Accordingly, the following sticking probabilities are used: P a i . j where the 
superscript “<z” indicates adsorption (to distinguish the adsorption process from the other processes), the left subscript index 
represents the flux species (denoted specifically by the indices s, hc, ao), and the right subscript index represents the 
surface species (denoted specifically by the indices g, s, hc, ao). Using these sticking probabilities, the fractional coverage of 
a species “f ’ on species is given by 


G°-j=l~e 


In) 


[5] 
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Fig. 1 - Two-level fractional coverage of unit area. 


Fig. 2 - Fractional coverage versus fluence for various P s . 


The quantity G represents a percentage of surface species j that is converted to species / and, therefore, is a relative fraction 
of Fj . This representation holds for the other processes to be discussed. Using adsorption as the example, it is emphasized 
that the relative coverage produced is due to a fluence m of species type i striking type j having n sites per unit area; the site 
density n includes all sites and remains fixed regardless of specie occupancy. Consequently, the resulting relative fractional 
coverage represents a change in absolute fractional coverage, brought about by the process occurring during the time interval 
associated with the prescribed fluence increment. The relative coverage brought about by the other processes is discussed in 
the following subsections. 

4.2 Desorption 

Since desorption is an activated process, it is dependent on substrate temperature. It is also dependent on surface coverage. 
As a substrate is heated, the rate of desorption may have a maximum due to greater difficulty in removing the more tightly 
bound molecules remaining at lower coverage. However, for a fixed coverage, the rate of desorption will generally be an 
increasing function of temperature. In the modeling, the processes are calculated assuming a fixed temperature. The fraction 
of type i molecules being desorbed from type j molecules in a unit time is denoted by G*_y. One model used here for 
desorption is given by 


Gflj = e~ cAhdllRT , [ 6 ] 

where c is a constant depending on the adsorbate and adsorbent combination. In this model, desorption occurs from the upper 
layer only. 

4.3 Reaction depletion by collision 

Atoms and molecules can also be removed from the surface due to reaction depletion. For example, AO striking hydrocarbon 
sites can cause the hydrocarbon to be activated sufficiently to escape the surface. This depletion mechanism is assumed to 
hold for all combinations except fluxes hitting glass sites, which are assumed inert to removal. The probability of type j 
molecules being depleted by collisions with type i flux in a unit time is denoted by G rc Vy. Since the statistics for depletion 
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follow identically from the statistics of adsorption, the expression for the depletion fraction is similar to equation [5] and is 
given by 


G?-j=l~e 


In) 


[7] 


where P n \.j is the probability that a site is depleted given an /-type fluence m, and overall site density n. Also the superscript 
“rJ” indicates reaction depletion by surface flux. Again, the left subscript index i represents the flux species (z = ao , s , and 
he ), and the right subscript index j represents the active surface species (j = ao , s , he). 

4.4 Chemical reaction 


For simplicity of modeling, only chemical reaction between AO and silicone are considered. Thus the remaining processes 
affecting surface coverage are: (1) chemical reactions due to collisions with the AO and silicone surface constituents by the 
opposite member’s flux and (2) chemical reactions due to surface reactions between these two species at adjacent sites. The 
flux-on-surface reactions are treated similarly to the adsorption and depletion processes given by equations [5] and [7], 
respectively. Hence, the fraction of sites converted in a unit time by reactions of type (1) can be written as 


Qtg = Y_ e - (m ‘ P H ,n) 


[ 8 ] 


The superscript fck rg” indicates reactions forming glass, and the subscript indices denote the species AO or silicone, with the 
flux component being the left index. The type (2) chemical reactions are taken to be proportional to the product of the surface 
concentrations. This proposed model seems reasonable. A more general approach would be to use the product of the 
concentrations raised to powers that are characteristic of the reaction and consistent with the experimental behavior; the 
exponents of the concentrations may even be fractional. The actual reaction mechanism may be complex, involving a number 
of elementary reactions and reaction intermediaries. Lacking definitive values of the exponents for the overall reaction 
process, it is assumed that the fraction of sites undergoing surface reactions of [AO] + [S] —> [G] depends bilinearly on the 
surface concentrations of these two species. Hence, the surface-reaction depletion fraction, denoted by G rg is a non- 
dimensional constant analogous to the rate constant: 


G rg = k . [9] 

This states that a fixed fraction of the reaction species will be converted in a unit time starting from specified concentrations. 
Extension of the model can be made to include chemical reactions between all species types, such as AO reactions with 
hydrocarbons, by including the appropriate terms to account for those chemical processes. 

What has been done here is to focus on only the sites occupied by reactant types, and thus the fractional gains above are 
relative to the sites of the respective species. To obtain the corresponding fractional specie coverage with respect to the total 
aggregate, the various G-quantities must be multiplied by their fractional representations. This is done in the next section 
using the appropriate weights through the F,- and Fy. 


5. ITERATION EQUATIONS FOR SURFACE COVERAGE 

The iteration equations updating the absolute fractional coverage from the previous unit time step are now presented. This is 
merely a matter of scaling the relative fractions for each process, done by multiplying by the appropriate fractional weights 
using the F,- and Fy. 

The iteration equation for the fractional glass surface coverage, F g , is explained for each group of terms according to process, 
which can be readily identified by the superscript indices. The superscripts k and k-1 are used to denote the new and previous 
time levels, respectively. The remaining species equations follow similarly and no further explanation is given. Before 
detailing the process, it is noted that for each equation every term contains a G factor except the leading term F*‘yX 1. The 
occurrence of a term in an equation is balanced by the corresponding negative term in another equation. It results, as 
expected, that the following conservation of fractions holds between time steps: 
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[ 10 ] 


F k + F k + F k + F k = F k ~ l + F k ~ 1 + F k ~ l + F k ~ l - 1 

r g ^ r s ^ r hc ^ r ao r g + r s + r hc + r ao 1 • 

The iteration equations for the glass, silicone, hydrocarbon and atomic oxygen fractions are given by the following equations: 

f,‘ = -g;,_, -g£_,) + F I ‘- , r-;(G- + g£.,+g^,) 

+ T?k-\ r'hc rd i/^rd \ , 77 A'-l 17 ao / rd , /^rd , s~i rd \ 

t hc t g ( Lr s-hc + U hc-hc + U ao-hc ) + r a o r g ^s-ao +(j hc-ao +(j ac-ao) 


+ F k ~ l F s G ds + F k ~'F hc G ds + F k ~'F ao G ds + F k ~ l G rg + F k ' l G rg +? F k ~ l F k ~'G rg rm 

r g^ J s ^ r hc r g KJ hc^ r ao r g kJ ao^ r s kJ ao-s^ r ao ^ s-ao ^ s r ao U ’ L U 1 


F s k = F k ~ l (l - G r a l s - G a hc _ s - Gl _, )+ F k g -'Gl g + F h k c -'G:_ hc + F^G° s _ ao + F k ~ l (F/ - 1)(G^ + G rd c _ s + G«_ s ) 
+ F s hc iQtbc + g £* + c: d 0 _ hc ) + F a k ;'Fr (g;1 + G2_ + G ld _ ao ) 


+ F?-\f; -1)G? +F^- l F s hc Gt + F“F."GZ - F k ~ l F k ~ l G rg , 


he ' * ao 


[ 12 ] 


K = K~ l (l - G a s . hc - G ° ao _ hc )+ F k -'G a hc _ g + F k -'G a hc _ s + F k 0 - l G a hc _ a0 + F^F,* (G* + G rd _ s + G r a d 0 _ s ) 
+ F?(F* - 1 )(G;^ + G£* + G£* ) + F fl r F* (G^ + G A t fl0 + G ' d _ ao ) 


+ F/-‘F; c Gf + F k ~ l (F h h c c - \)Gt + F“FZG£ , 


he ao he ^ ao 


[13] 


f.‘ = C (i - G2.„ - G‘,-„ - c;U )+ + FFGI-, + C G L-* + f;-' f; ( c;', + g£, + at , ) 

+ C C «£* + G* + G* ) + F*- 1 (F" - 1)(G" + at _ + G- ) 


. /r^- 1 Z7 5 # 7 * + r k ~ { r hc ri ds 4 . F k ~ l (F ao — 1 

^ r s r ao KJ s ^ r he r ao KJ hc^ r ao \ r ao l > KJ an r .v r ao U 


[14] 


In discussing the terms in equation [1 1], the terms are examined from left to right. The first terms involve the adsorption of 
silicone, hydrocarbon and AO onto glass. The increment to absolute fraction of glass by adsorption is negative due to loss of 
glass fraction. The increment equals the product resulting from multiplying the sum of the three G a i. g terms by the absolute 
fractional area of glass, F*' 1 g existing at the previous time step. The second category comprises three groups of three terms. 
These consist of the reaction depletion caused by atomic and molecular flux dislodging species from their sites. Since it is 
assumed that silica sites are not removable, the terms are positive here due to this process. The first set of three depletion 
terms is the set due to silicone being removed by silicone, hydrocarbon, and AO flux, respectively. The absolute fraction of 
glass surface gained is dependent proportionally on the absolute fraction of silicone, times the relative fraction of this 
silicone fraction that covers glass, P g . Thus, the product of these factors, multiplied by the sum of the three depletion terms 
in the silicone group, gives the increment due to silicone depletion. The two remaining depletion groups involve the 
hydrocarbon and AO site fractions, and they follow similarly by using the appropriate multiplicative fractions. The third 
group consists of three terms due to desorption. The first of these involves desorption of silicone. The absolute fractional gain 
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in glass area due to silicone desorption is equal to the triple product of the absolute fraction of silicone, F*‘ y s , times the 
relative fraction of silicone that covers glass, times the fraction of silicone desorbed in the time period, G ds s . The two 
desorption terms relating to hydrocarbon and AO follow similarly. The fourth set of terms represent the creation of glass area 
by silicone and AO flux and are positive terms. These involve the G rg ao . s and G rs s . ao relative fractions, which must be 
multiplied, respectively, by F*' y 5 and F*' 1 ao representing the absolute fractions for the surfaces struck by the opposite type 
flux. The last group consists of the singleton term due to surface reaction of silicone and AO. It is equal to the product of the 
absolute silicone, F 4 "' 1 s , and AO, F*' 1 ao , fractions times the twice the fraction, G' ;? , of either of these species that is converted in 
a time step. 


6. ITERATION EQUATIONS FOR SUB-LEVEL COVERAGE 

The iteration equations (with k denoting the time-step level) for the glass, silicone, hydrocarbon and atomic oxygen absolute 
fractions,/^, are given by the following equations: 

6.1 Glass on glass 


g k e k-1 k-1 g k-1 a a a k-1 s k-1 g k-1 rd rd rd ds rg 

(f)=(f) -F (F) [G + G + G ] + F (F) {(F) [G + G +G +G ] + G 

v g 7 v g 7 g g 7 L s-g hc-g ao-g J s x g 7 1 g s-s hc-s ao-s s ao-s 7 

k-1 he k-1 g k-1 rd rd rd ds k-1 ao k-1 g k-1 rd rd rd ds rg 

+ F (F ) (Fj [G . + G . . + G + G . ] +F (F ) (F ) [G +G h +G +G +G 

he g g s-hc hc-hc ao-hc he ao g g s-ao hc-ao ao-ao ao s-ao 

k-1 k-1 s ao rg 

+ F F [F + F ]G [15] 

s ao g g 


6.2 Glass on silicone 


g k g k-1 k-1 g k-1 a a a k-1 s k-1 g k-1 rd rd rd ds s k-1 rg 

(f ) = (r ) - F (F ) [G + G + G ] + F {(F ) (F ) [G + G . +G + G ] + (F ) G } 

v s 7 v s 7 g s 7 s-g hc-g ao-g s 1 g s L s-s hc-s ao-s s s ao-s 


g s' s-g hc-g 

k-1 he k-1 g k-1 rd rd 


ao-g 
rd 


g' s' s-s hc-s ao-s 

,.w « . b ds " k-1 ao k-1 g k-1 rd rd rd ds 

+ F u (F ) (F ) [G + G + G +G1+F {(F ) (F ) [G + G +G +GJ 

he 8 s s-hc hc-hc ao-hc he ao g s s-ao hc-ao ao-ao ao 

ao k-1 rg k-1 k-1 s ao rg w P4 „ 

+ (F ) G }+ F F [F + F ]G [16] 

s s-ao 7 s ao s s 


6.3 Glass on hydrocarbon 

g k g k-1 k-1 g k-1 a a a k-1 s k-1 g k-1 rd rd rd ds s k-1 rg 

(f ) = (f ) - F (F ) [G + G + G ] + F {(F ) (F J [G + G +G + G ] + (F ) G 

v he 7 v he 7 g he L s-g hc-g ao-g J s 1 g 7 he s-s hc-s ao-s s he ao-s 

k-1 hc k-l g k-1 rd rd " rd ds k-1 ao k-1 g k-1 rd rd rd ds 

+ F (F ) (F\ ) [G + G +G +G ] + F {(F ) (F ) [G +G ] 

v O 7 v hr 7 L c-hr hr-hr an-hc hc J ao 1 e he S-ao hc-ao ao-ao ao 

[17] 


s-hc hc-hc ao-hc he 

ao k-1 rg k-1 k-1 s ao rg 

+ (F ) G }+ F F [F +F 1G 

h '' s-ao s ao 


6.4 Glass on AO 

k ^ k-1 _ k-1 g k-1 


g k g k-1 k-1 g k-1 a a a k-1 s k-1 g k-1 rd rd rd ds s k-1 rg 

(f ) = (f ) - F (F ) [G +G. +G ] + F {(F) (F ) [G + G +G + G ] + (F ) G } + 

v ao 7 v ao 7 g v ao 7 L s-g hc-g ao-g s 1 g ao s-s hc-s ao-s s ao . ao * s 

k-1 hc k-l g k-1 rd rd rd ds k-1 ao k-1 g k-1 rd rd rd ds 

F (F)(F)[G + G + G + G ] + F {(F)(F)[G ] 

he v g 7 v ao 7 L s-hc hc-hc ao-hc he ao g ao s-ao hc-ao ao-ao ao 

ao k-1 rg k-1 k-1 s ao rg 

+ (F ) G~ }+ F F [F + F ]G [18] 

ao s-ao 7 s ao ao ao 

6.5 Silicone on glass 


s k s k-1 k-1 a k-1 s k-1 s k-1 rd rd rd ds a a a rg 

(f ) =(f ) + F G -F (F ) { ( 1-(F ) ) [G + G +G +G] + G +G h +G +G } 

v z v g 7 g s-g s v g 7 1 v v s 7 7 L s-s hc-s ao-s s s-s hc-s ao-s ao-s 

k-1 he k-1 s k-1 rd rd rd ds k-I ao k-1 s k-1 rd rd rd ds 

+ F (F ) (F ) [G +G..+G +G1 + F (F)(F)G + G + G +G] 

he v s 7 v g 7 L s-hc hc-hc ao-hc hc J ao v s g s-ao hc-ao 9A - an 

k-1 k-1 s rg 

-F F F G 

s ao g 


[ 19 ] 
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6.6 Silicone on silicone 


s k s k-1 k-1 s k-1 a s k-1 s k-1 rd rd rd ds a a re 

(f)=(f) +F { (1-(F ) ) G -(F) [( 1-(F ) ) [G + G + G +G ] + G +G +G* ]} 

s s s s s-s s s s-s hc-s ao-s s hc-s ao-s ao-s J 

k-1 he k-1 s k-1 rd rd rd ds k-1 ao k-1 s k-1 rd rd rd ds 

+ F (F) (F) [G +G +G + G ] + F (F) (F) [G + G , + G + G ] 

he s s s-hc hc-hc ao-hc he ao s s s-ao hc-ao ao-ao ao 

k-1 k-1 s rg 

- F F F G‘ [20] 

C Jin C L J 


s-hc 


6.7 Silicone on hydrocarbon 


s k s k-1 k-1 s k-1 s k-1 rd rd rd ds a a a rg 

( f J=(f J -F <FJ 1KHFJ ) [G „ + G + G + G ] + G + G +G + G } 


he he 7 s he 1 v s' ' L s-s hc-s ao-s s J ' s-s hc-s ao-s ao-s' 

k-1 a he k-1 s k-1 rd rd rd ds k-1 ao k-1 s k-1 rd rd rd ds 

+ F [G + (F ) (F ) [G + G + G + G ]] + F (F ) (F ) [G + G + G + G ] 

he s-hc s he s-hc hc-hc ao-hc he ao s he s-ao hc-ao ao-ao ao 

k-1 k-1 s rg 

-F F F G [21] 

s ao he J 


s-hc 


6.8 Silicone on AO 


s k s k-1 k-1 s k-1 s k-1 rd rd rd ds a a a rg 

(f )=(f ) -F (F ) [(1-(F ) ) [G + G + G + G ] + G +G +G +G ] 

ao ao s ao s s-s hc-s ao-s s s-s hc-s ao-s ao-s 

k-1 he k-1 s k-1 rd rd rd ds k-1 a ao k-1 s k-1 rd rd rd ds 

+ F (F ) (F ) [G k + G + G +G1 + F [G + (F ) (F ) [G + G + G + G ]] 

" u " u ~ ao-hc he ao s-ao s ao s-ao hc-ao ao-ao ao 

k-1 k-1 s rg 

-F F F G [22] 

s ao ao 


s-hc 


hc-hc 


6.9 Hydrocarbon on glass 


he k he k-1 k-1 a k-1 s k-1 he k-1 rd rd rd ds 

(f ) = (f ) + F G + F (F ) (F ) [G + G + G + G ] 

g g g hc-g s he g s-s hc-s ao-s s 

k-1 he k-1 he k-1 rd rd ra ds a a a 

- F (F ) [ (1-(F ) ) [G + G + G +G1 + G +G +G J 

he g he s-hc hc-hc ao-hc he s-hc hc-hc ao-hc 

k-1 ao k-1 he k-1 rd rd rd ds 

+ F (F fc ) (F ) [G + G . + G + G ] 

ao he g s-ao hc-ao ao-ao ao 


[23] 


6.10 Hydrocarbon on silicone 


he k he k-1 k-1 a s k-1 he k-1 rd rd rd ds 

(f)=(f) + F [G + (F ) (F ) [G +G.+G +G1] 


hc-s he 7 s 7 s-s hc-s ao-s 

k-1 he k-1 he k-1 rd rd rd ds a 

F. (F ) [ (1-(F ) )[G + G + G +G1 + G +G +G J 

s he s-hc hc-hc ao-hc he s-hc hc-hc 

k-1 ao k-1 he k-1 rd rd rd ds 

+ F (F J (F J [G . „ + G k _ „ + G _ + G ] 


he 


ao-hc 


he 


hc-ao 


[24] 


6.11 Hydrocarbon on hydrocarbon 


he k he k-1 k-1 s k-1 he k-1 rd rd rd ds 

( f J = < f J + F <F „ ) ( F k > [° + G . + G + G ] 

he he s he he s-s hc-s ao-s s 

k-1 he k-1 a he k-1 he k-1 rd rd rd ds 

+ F (1-(F ) )G -(F ) [(1-(F . ) ) [G + G + G +GJ 

he he hc-hc he he s-hc hc-hc ao-hc hc J 

a a k-1 ao k-1 he k-1 rd rd rd ds 

+ G +G ]} + F (F ) (F ) [G + G + G +G] 

s-hc ao-hc ao he he s-ao hc-ao ao-ao ao 


[25] 


6.12 Hydrocarbon on AO 


he k he k-1 k-1 s k-1 he k-1 rd rd 

(f ) =(f ) +F (FJ (F ) [G + G 


rd ds 

x , v- XL- . — . + G +G] 

uv/ uv/ s he 30 s-s hc-s ao-s s 

k-1 he k-1 he k-1 rd rd rd ds a a a 

- F (F ) [(1- (F ) ) [G + G + G +G J + G + G + G J 

he ao he s-hc hc-hc ao-hc he s-hc hc-hc ao-hc 

k-1 a ao k-1 he k-1 rd rd rd ds 

+ F [G + (F ) (F ) [G + G + G + G ]] 

ao hc-ao he ao s-ao hc-ao ao-ao ao 


[26] 
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6.13 AO on glass 


ao k ao k-1 k-1 a k-1 s k-1 ao k-1 rd rd rd ds 

(f ) = (f ) + F G + F (F ) (F ) [G + G + G + G ] 

e e g ao-g s ao g s-s hc-s ao-s s 

k-1 he k-1 ao k-f rd w rd rd ds k-1 ao k-1 ao k-1 rd rd rd ds 

+ F (F ) (F ) [G +G +G + G J- F (F ) [(1-(F ) ) [G +G k +G + G ] 

he ao g s-hc hc-hc ao-hc he ao g ao s-ao hc-ao ao-ao ao 

a a a rg k-f k-1 ao rg 

+ G + G + G + G ] - F F FG [27] 

s-hc hc-ao ao-ao s-ao s ao g 

6.14 AO on silicone 


ao k ao k-1 k-1 a s k-1 ao k-1 rd rd rd ds 

(f ) = (f ) + F [G + (F ) (F ) [G + G + G + G ] 

s s s ao-s ao s s-s hc-s ao-s s 

k-1 he k-1 ao k-1 rd rd rd ds k-1 ao k-1 ao k-1 rd rd rd ds 

+ F (F ) (F ) [G +G +G + G ] - F (F ) [(1- (F ) ) [G + G +G +G 

he ao s s-hc hc-hc ao-hc he ao s ao ark 

a a a rg k-1 k-1 ao rg 

+ G + G + G + G ] - F F F G 

s-ao hc-ao ao-ao s-ao s ao s 


s-ao hc-ao ao-ao ao 


[28] 


6.15 AO on Hydrocarbon 

ao k ao k-1 k-1 s k-1 ao k-1 rd rd rd ds 

(f ) = (f ) + F (F ) (F ) [G + G + G + G ] 

he he s ao he s-s hc-s ao-s s 

k-1 a he k- 1 ao k- 1 rd rd rd ds 

+ F_ [G + (F ) (F ) [G +G +G + G ]] 

he ao-hc ao he s-hc hc-hc ao-hc he 

k-1 ao k-1 ao k-1 rd rd rd ds a a a rg k-1 k-1 ao rg 

-F (F . ) [(1- (F ) ) [G + G + G + G ] + G +G k +G +G -F F FG [29] 

ao he ao s-ao hc-ao ao-ao ao s-ao hc-ao ao-ao s-ao s ao he 

6.16 AO on AO 

ao k ao k-1 k-1 s k-1 ao k-1 rd rd rd ds 

(f ) =(f ) +F (F ) (F ) [G + G + G + G ] 

ao ao s ao ao s-s hc-s ao-s s 

k-1 he k-1 ao k-1 rd rd rd ds k-1 ao k-1 a ao k-1 ao k-1 

+ F. (F ) (F ) [G + G + G +G 1 + F (1-(F ) )G -(F) (1-(F ) ) 

he v ao y v ao L s-hc hc-hc ao-hc he ao ao ao-ao ao ao 

rd rd rd ds a a rg k-1 k-1 ao rg 

[G + G + G + G ] + G + G + G ] } - F F FG [30] 

s-ao hc-ao ao-ao ao s-ao hc-ao s-ao s ao ao 


These updated absolute fractions, are used in computing new relative fractions. F j, by using Equation 3. 


7. RESULTS OF MODELING 

For the purpose of examining the essential features of film growth, a prototype code was constructed using a popular 
spreadsheet computer program to simulate the processes of adsorption, desorption, depletion, and chemical reaction. Using 
this code, SiO x growth can be modeled based on the relative fluxes of AO, hydrocarbon, and silicone. By knowing the 
amount of contaminants trapped during the film growth, the resulting optical transmittance of the contaminated film can be 
modeled based on assumed values of absorption coefficients. The preliminary results reported here are based on simple 
models that are hoped to be useful in explaining results from in-space and ground-based contamination. A brief review of 
some recent findings concerning both in-space and ground-based contamination is presented before discussing results from 
this present study. 

In-space contamination results 5 indicate that the flux of volatile components of materials such as silicones and hydrocarbons 
decreases over time from the early days of a mission. From studying the variation of thickness among the different layers in 
films deposited on spacecraft material exposed to both contaminant and AO flux, it was concluded that variations in the rate 
of oxidation occurred. The thickness of the layers decreased from the inner layers to the outer layers of the films. Variations 
in the degree of oxidation were also determined by examining surfaces exposed to different rates of contaminant flux. The 
oxidation was determined to vary with respect to the rate of silicone and hydrocarbon arrival because a much thinner, more 
fully oxidized and more transparent contaminant layer was formed on the surfaces believed to have received lower silicone 
and hydrocarbon flux. Ground-based results 5 support the conclusion that small fluxes of contaminants relative to a fixed AO 
flux produce a nearly fully oxidized film, whereas large ratios of contaminant to AO flux result in markedly reduced 
oxidization fractions. 
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A test case is presented using the numerical code and is compared to the qualitative results observed experimentally. The 
specification of a decaying contaminant flux is made because, as mentioned above, the flux of volatile materials diminishes 
with time. An AO flux composed of periodic spikes of large amplitude superimposed on a relatively low base level is 
prescribed. First, the variation in fractional coverage by species is discussed, and then a simple model of optical transmittance 
degradation is presented along with a result illustrating the effect of the AO-to-silicone flux ratio. 

The distributions of the silicone and hydrocarbon fluxes for the test are shown in figure 3 as solid curves. The AO flux, 
having a peak relative amplitude of one, is shown as the dotted curve. Based on these fluxes, a somewhat arbitrary or 
hypothetical specification of values for the various physical quantities such as sticking probability was made. The actual 
values used are values chosen as being representative in the sense that they are reasonable, but not necessarily and probably 
not correct. For example, a sticking probability equal to 0. 1 may be reasonable for a single gas-solid specification but not 
likely the correct choice for all species, substrate temperatures, and fractional coverage. For simplicity in specification, 
essentially uniform probabilities were used for the various processes. In any case, the various probabilities must be selected 
in conjunction with the prescribed fluxes so as to limit the change in fractional surface areas during a unit time step. Although 
not mentioned previously, the procedure relies on the assumption that the processes can be treated independently. This will 
not be the case if the average site is being multiply struck in a single time-step increment. This restriction is easily met by 
using relatively small m-lo-n ratios. For such a set of hypothetical probabilities, figure 4 shows the variation in surface 
coverage by species corresponding to the relative fluxes presented previously. The top curve represents the fractional 
coverage of glass starting from a clean silica surface. From this figure, it can be seen that the glass fractional area initially 
declines quickly as the other species populate through adsorption. The hydrocarbon and silicone fluxes increase steadily at 
first, peak and then begin to decline in coverage. This decline in contaminant coverage is balanced by an increase in glass 
coverage and a slight increase in AO coverage. The jumps in the AO coverage, obviously corresponding to the spikes in the 
AO flux, are balanced by almost equally sharp falls in the hydrocarbon coverage and silicone coverage, with slight drops in 
the glass coverage. Variation of parameters about values used in this example produce reasonable changes in the fractional 
coverage by species. 

One important component contributing to the fractional glass coverage is that which accrues due to the chemical reaction 
processes associated with the final three terms of Equation [11]. These three terms are the source for film growth. The first 
two terms represent the fractional increase in glass by chemical reactions resulting from flux of type i on surface type j\ 
where the species type is either silicone or AO. The third of these three terms represents the fractional increase in glass due to 
chemical reactions occurring at adjacent silicone and AO surface sites. A simple modeling of the contaminant film would be 
to use a linear combination of the three terms while accounting for trapping of the contaminants. To simplify the discussion, 
the choice of using only the last term of the three is made here and, furthermore, the assumption is made that trapping of 
contaminants is proportional to this term. Using the fractional coverage provided in the code, the contaminant thickness can 
be integrated numerically. The relative transmittance is calculated based on the Beer-Lambert law, which is given, assuming 
the absorbing material is distributed uniformly throughout the film, by 

r = <r°*\ [31] 

where a is the optical absorption coefficient and L is the integrated film thickness, which generally increases nonlinearly 
with time. The time variation of relative transmittance of a contaminant film produced using the fluxes as prescribed above is 
presented in figure 5 as the solid curve. The rate of reduction in transmittance increases significantly during intervals of peak 
AO flux. The exponential-like decline in transmittance is due to the exponentially declining silicone flux. The variation of 
transmittance is also shown for the same case but where the AO flux is prescribed as being a constant, equal to the time- 
averaged value of the spiked AO-flux distribution. The transmittance decline seen for the constant AO-flux case is slightly 
less than the spiked AO-flux case, otherwise it nearly serves as an upper envelope for the spiked AO-flux transmittance 
curve. The greater overall decline in transmittance for the spiked AO-flux distribution can be explained by noting that the 
early large spikes in AO flux provide the oxidizer needed to combine with the high initial contaminant flux. This produces a 
rapid decline in transmittance that is never again matched. For “rear values of the parameters, this would indicate that the 
most significant contaminant degradation may occur relatively early in a spacecraft mission. It must be remembered that the 
values used here for the probabilities of adsorption and reaction, etc., are reasonable hypothetical values and are mainly used 
as placeholders to check out the mechanics and the features of the code. The correct values of the process parameters must 
come from either experiment or thermodynamic analysis. 


NASA/TM — 200 1 -2 1 0597 


10 



Fig. 3 - Variation of species flux with time. 


Fig. 4 - Evolution of fractional coverage for species. 
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8. SUMMARY 


Contaminant films can form on spacecraft surfaces degrading system performance. AO present in LEO can react with volatile 
silicone and hydrocarbon molecules to produce a glassy SiO x based contaminant coating. In order to model the growth of the 
contaminant film, procedures were devised to model the surface processes such as adsorption, desorption, collision depletion, 
and chemical reaction. It was found that a probabilistic approach could be used to estimate the likelihood of a process in 
terms of specified fluxes and existing surface coverage. Four specie types, consisting of silica glass, silicones, hydrocarbons, 
and AO, were considered. Because the different processes involved changes in surface site occupancy, it proved convenient 
to conceptually group sites by specie type. A detailed bookkeeping method devised to track surface coverage by specie type 
allowed the surface fractional coverage and sub-layer fractional coverage to be updated based on the changes produced 
during fluence increments. Equations were presented for updating the fractional species coverage based on changes in specie 
site occupation, including a full discussion for one such equation. Also, equations were presented, according to each 
overlaying specie type, for the first sub-layer fractional coverage by specie. Because the values of the various parameters 
characteristic of a given set of flux distributions are not generally known, estimated values must be used in many cases. A 
standard spreadsheet program proved a suitable platform from which to construct a prototype numerical code. Example 
results using a code that calculates the evolution of surface coverage were presented. One example illustrates for a particular 
case how the fractional coverage evolves given time dependent fluxes. An example for a simple modeling of transmittance 
indicates that the growth in a contaminant layer may come early in a spacecraft mission if the volatile contaminant fluxes 
diminish exponentially with time. 


9. REFERENCES 

1. B.A. Banks, J.A. Dever, L. Gebauer, and C.M. Hill: “Atomic Oxygen Interactions with FEP Teflon and Silicones on 
LDEF\ Proceedings of the 1 st LDEF Post-Retrieval Symposium , pp. 801-816, Kissimmee, Florida, June 2-8, 1991. 

2. E.R. Crutcher and K.J. Warner, “Molecular Films associated with LDEF’, Proceedings of the 1 st LDEF Post-Retrieval 
Symposium , pp. 155-177, Kissimmee, Florida, June 2-8, 1991. 

3. Y. Noter, E. Grossman, L. Genkin, Y. Haruvey, M. Murat, I. Ross, R. Vered, Y. Lifshitz, Y. Efraty, H. Granot, 
A. Halsbersberg, H. Lauria, and A. Stoler, “Variations in the Telemetry of an Offeq Satellite’s Sun Sensors”, 
Proceedings of the 3 rd International Space Forum on Protection of Materials and Structures from the LEO Space 
Environment , Toronto, Canada, April 25-26, 1996. 

4. B. Blakkolb, L.E. Ryan, H.S. Bowen and T.J. Kosic: “Optical Characterization of LDEF Contaminants”, Proceedings of 
the 2 nd LDEF Post-Retrieval Symposium , pp. 1035-1040, San Diego, California, June 1-5, 1992. 

5. B.A. Banks, S. Rutledge, E. Sechkar, T. Stueber, A. Snyder, K.K. de Groh, C. Haytas, and D. Brinker: “Issues and 
Effects of Atomic Oxygen Interactions with Silicone Contaminant on Spacecraft in Low Earth Orbit”, Proceedings of the 
8 th International Symposium on Materials in a Space Environment , Arcachon, France, June 5-9, 1991 

6. G. Attard and C. Barnes, Surfaces , Chapter 1, Oxford University Press, Oxford, 1998. 

7. A.W. Adamson, Physical Chemistry of Surfaces, Wiley Interscience, New York, 1990. 


N AS A/TM— 200 1 -2 1 0597 


12 


REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY ( Leave blank) 

2. REPORT DATE 

June 2001 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 

5. FUNDING NUMBERS 


Modeling of Transmittance Degradation Caused by Optical Surface 
Contamination by Atomic Oxygen Reaction With Adsorbed Silicones 


6. AUTHOR(S) 

Aaron Snyder, Bruce Banks, Sharon Miller, Thomas Stueber, and Edward Sechkar 


WU-755-A4-1 3-00 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E- 12555 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM— 2001-210597 


11. SUPPLEMENTARY NOTES 

Prepared for the 45th Annual Meeting, the International Symposium on Optical Science and Technology sponsored by the 
International Society for Optical Engineering, San Diego, California, July 30-August 4, 2000. Aaron Sfciyder, 

Bruce Banks, and Sharon Miller, NASA Glenn Research Center; Thomas Stueber and Edward Sechkar, Dynacs 
Engineering Company, Inc., 2001 Aerospace Parkway, Brook Park, Ohio 44142. Responsible person, Aaron Snyder, 
organization code 5480, 216-433-5918. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category: 18 

Available electronically at http://gltrs.grc.nasa.gov/GLTRS 

This publication is available from the NASA Center for AeroSpace Information. 301-621-0390. 


Distribution: Nonstandard 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

A numerical procedure is presented to calculate transmittance degradation caused by contaminant films on spacecraft 
surfaces produced through the interaction of orbital atomic oxygen (AO) with volatile silicones and hydrocarbons from 
spacecraft components. In the model, contaminant accretion is dependent on the adsorption of species, depletion reactions 
due to gas-surface collisions, desorption, and surface reactions between AO and silicone producing SiO x (where x is 
near 2). A detailed description of the procedure used to calculate the constituents of the contaminant layer is presented, 
including the equations that govern the evolution of fractional coverage by specie type. As an illustrative example of film 
growth, calculation results using a prototype code that calculates the evolution of surface coverage by specie type is 
presented and discussed. An example of the transmittance degradation caused by surface interaction of AO with deposited 
contaminant is presented for the case of exponentially decaying contaminant flux. These examples are performed using 
hypothetical values for the process parameters. 


14. SUBJECT TERMS 

Contamination; Low Earth orbit; Atomic oxygen; Silicones; Sensor; 
Adsorption; Optical transmittance 


15. NUMBER OF PAGES 

18 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 


